Electron transmission in normal/heavy- fermion superconductor junctions 
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The Andreev reflection between a normal metal (N) and a heavy-fermion superconductor (HFS) 
is studied and the boundary conditions for the electron's wave function in the two systems are 
established in the framework of a two band model for the HFS. Hence we show in a simple and 
explicit way that the mass enhancement factors in the heavy-fermion (HF) metal do not cause 
impedance at the N/HFS interface, in accordance with arguments previously presented. We also 
present an extension of the theory to a two-fluid model for the heavy-fermion, as possibly applicable 
to e.g., CeCoIns. 
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I. INTRODUCTION AND MODEL 

Electronic scattering at the interface between a nor- 
mal metal (N) and a superconductor has been studied 
by Blonder, Tinkham and Klapwijk 1 (BTK) and it was 
found that the subgap conductance is enhanced due to 
Andreev reflection. On the other hand, the Fermi veloc- 
ity mismatch between the two metals always produces 
an effective barrier which decreases the conductance^. 
If the superconductor is a heavy-fermion (HFS), the 
greater Fermi velocity mismatch would lead us to expect 
a strongly reduced subgap conductance. Experimentally, 
however, the subgap conductance does not seem to be 
strongly reduced in N/HFS junctions. An argument to 
explain this behavior has been put forward by Deutchcr 
and Nozieres 3 , who claimed that mass enhancement fac- 
tors in the heavy-fermion metal do not cause impedance 
at the interface. 

Motivated by recent experiments^ on Au/CeCoIns in- 
terfaces, we study electron scattering at the interface be- 
tween a normal (light) metal (N) and a heavy-fermion 
superconductor. Starting from a more realistic two-band 
model for the HFS, where a conduction c-electron band 
hybridizes with a localized /-electron band, plane-wave 
solutions for each bulk subsystem can be written down. 
We explicitly obtain the matching conditions for the wave 
function at the interface and confirm the claim in Rcf 3 , 
by explicitly showing that the electron velocities involved 
are those of the conduction c-bands. Since the localized 
/-electrons are dispersionless, only the c-conduction elec- 
trons satisfy matching conditions at the interface. This 
is the basis of the result. 

An extension to a two-fluid model of the HFS, where a 
normal light fluid coexists with the superconducting one, 
is also presented, which may be relevant to e.g., CeCoIns. 
Our theory is analogous to the quantum waveguide the- 
ory for mesoscopic structures^. 

The model for the HFS is: 
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The superconducting pairing, A k , has been explicitly 
written for the /-electrons because such a model qualita- 
tively reproduces the generic phase diagram of a HFS^^. 
Equation fl} is an effective model obtained in the large 
Anderson-?/ limit, where it is understood that both the 
/-electron level, et, and the hybridization, V, are ob- 
tained self-consistently. 

The Bogolubov operators which diagonalize the above 
Hamiltonian are given by: 



7kcr = uckcr - crv*ff_ k _ a + uf k(7 - cru*/i k _ CT 
where the amplitudes obey the linear system 
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The excitation energies are given by: 
25(k) 
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2 7 (k) = ( £/ - A ,) 2 + ( ek - M ) 2 + |A k | 2 + 2y 2 , 
/3 2 (k) = (e k - M ) 2 |A k | 2 +[(e / -M)(e k - M )-^ 2 ] : 
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For convenience, we introduce the simplified notation: 
£ = £k - /« , e = tf-fi, A = A k . (8) 
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The linear system ([3])-((4]) yields a relation between the 
amplitudes on the / and c subsystems: 
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From ^ we see that the amplitudes on the / sites are 
not independent being proportional to the amplitudes on 
the c sites. The description of the bands is standard and 
is briefly reviewed in the Appendix. 



II. ELECTRON TRANSMISSION 

An incident electron from the light metal, with energy 
E > measured from the Fermi level, will penetrate the 
HFS. The coherence factors for the transmitted quasi- 
particles (quasi-hole and quasi-electron) can be obtained 
from the linear system ©-dU). Using ^ to eliminate u, 
v we get: 

[V 2 - (E-£)(E-ej\ u- A- {E + £)v = 

A- (£-£>+ [V 2 - (E + £,)(E + e)] v = (10) 

Equating the corresponding determinant to zero we ob- 
tain: 



eV 2 ± V-A 2 !/ 4 + E 2 (e 2 + A 2 + V 2 - E 2 ) 
e 2 + A 2 - E 2 



(11) 



Equation (jlip determines the momenta of the transmit- 
ted quasiparticles (equation (Tl9|) below). 

At the Fermi level (E = 0), the quasi-particles decay 
exponentially into the HFS if the argument of the square 
root is negative (roughly if E < A*., for a specific direc- 
tion). Equation (fTTj) then gives 
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where we have defined 
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It is seen from equation (|13[) that the momentum on the 
decaying quasiparticle in the HFS has a real part which 
is kp and an imaginary part, k. The decay length, K , 
is determined by the slow Fermi velocity of the heavy- 
fermion system , v~p f 2 jV 2 . 

We write the electron wave function in the N/HFS sys- 
tem as a four-component vector by generalizing the col- 
umn vector in equation ([3]): 
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Inside the normal (light) metal we have / = g = 0, of 
course. We identify the N/HFS interface with the x = 
plane. When a plane wave is scattered at the interface, 
the wavevector component that is parallel to the interface 
is conserved. While no boundary conditions are imposed 
on / and g, the functions / and g are required to be con- 
tinuous at the NS boundary. Allowing for two different 
c-band effective masses in the two metals, we write the 
kinetic energy operator as&i£ 
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where m{x) = m n in the light metal (x < 0) and m(x) = 
rrihf in the heavy metal {x > 0). The mass rrihf is not 
large: it is simply the c-band effective mass. Because of 
the above form of the kinetic energy, the derivative of ip 
is not continuous at x = Q. The functions /, / and g, g 
are coupled by the equations: 



Ef = T-fi + US(x) f + VQ{x)f 



(15) 
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T + li-U5(x) g-VQ(x)~g, (16) 



where U denotes a potential barrier at the interface 
and O(x) denotes Heaviside's function. The remaining 
boundary condition for / and g is obtained upon inte- 
gration of (fl"5 |) -(fl"6 |) between x = CP and x = + : 
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where the prime denotes derivation with respect to x. 

We write ^(r) = e tk ^' v ^ip(x) where k|| and r 1 1 are the 
wavevector and position vector components parallel to 
the interface, and 



ip(x < 0) 
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where b denotes the particle-particle reflection amplitude 
and a denotes the Andreev reflection amplitude. The 
transmitted quasi-electron and quasi-hole in the HFS 
have amplitudes c and d, respectively. The momenta 
are obtained from energy conservation: E(k\\,k~) — 

E{k\ b k+) = -Uk\\,p-) = Uk\\,P + ) = E, where E{k) 
is given in equation (|A5|) and ^ n (k) denotes the incident 
electron dispersion in the normal metal measured from 
the Fermi level. 

In the following we calculate the reflection amplitudes 
a and b for E — and normal incidence (pu = 0), for 
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simplicity. In this case all functions depend on the co- 
ordinate x, only. Equations (fTT)| - (fTTj) give v + = —iu 
and k + = kp + in for the transmitted quasi-electron, and 
V- = iu and k~ = kp — in for the transmitted quasi-hole. 
Owing to the different c band effective masses, we denote 
the Fermi velocity in the light metal by v n , which should 
be comparable to that obtained from the slope of the c 



band in the HF, 



.,( c ) 



The continuity of the functions / 



and g at the interface implies: 

I + b = cu- 
ia = cu 

and equations (fTT)) - (fT5)) imply: 
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where fvpp denotes the Fermi momentum in the normal 
metal, v n — hpp/m n and vp^ = hkp /m^f. Introducing 
the dimensionless barrier parameter— 

z = f • 

nv n 

we may write the reflection amplitudes as: 
-2i ^~ 

(25) 
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The expressions ([25]) and (|26p are precisely what would 
be obtained if the HFS was a one-band superconductor 

(c) 

with a Fermi velocity v F . This can be traced back to 
the fact that only the c-electron parts of the wave func- 
tion, /(r), and g(r), satisfy matching conditions at the 
interface that are the same as in the case of a one-band 
superconductor. Inside the HFS, the /-site amplitudes 
/(r) and g(r) are directly proportional to the c-electron 
amplitudes as shown in @. The mass enhancement, or 
slow Fermi velocity, in the HF appears in (|2"5|) and (|2"6"1) 
through the coherence length n, as can be seen from (fT3"|). 
In the case of a clean interface (Z = 0) and long decay 
length (such as close to a nodal direction), k — > 0, these 
expressions simplify to 
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where r\ = Vp /v n . If there is no mismatch of light Fermi 
velocities (77 = 1 ) then 6 = and \a\ 2 = I leading in 
this rather special case to a perfect doubled conductance: 



f 



\b\ 2 
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III. TWO-FLUID HFS 

A theory for electron transmission from a normal metal 
to a one-band superconductor— was applied to an in- 
terface with the HFS CeCoIng^£. This material has a 
complex Fermi surface^ and seems to be well described 
by a two-fluid mode fe 13 i 14 . The superconducting gap 
has d-wave symmetry, as determined directly by various 
experiments ^ 15 ' 16 ! 17 . 

The HF metal CeCoIns is known to have at least one 
band of uncondensed light carriers in addition to the 
heavy fermion superconducting liquid. Recent Andreev 
reflection studies assume that the subgap conductance 
is a weighted average between that given by the BTK 
theory and a flat conductance due to the band of un- 
condensed fermions^. This motivates us to analyze the 
case where the incident electrons from the normal metal 
(labeled as system "I") can tunnel simultaneously to a 
c-band coupled to the / subsystem (labeled "2") and to 
a light uncondensed conduction band (labeled "3"), co- 
existing in the heavy-fermion material. Metal 1 is in the 
x < half-space and metals 2 and 3 are in the x > 
half-space, with the interface at x = 0. The theory we 
employ is an extension of the quantum waveguide theory 
of mesoscopic structures 6 - where a delta function poten- 
tial is introduced at the interface and where one of the 
circuit branches, "2", is the HFS model of equation ([TJ) . 

The potential U8(x — e) is in metal 1 {x < 0) and 
we shall take the limit e — > 0~. The wave function in 
the normal single-band metal 1 has particle and hole 
components^: 



ipi(x < e) 
and 
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The wave function in metal 2 is written in the same form 
as (HI: 
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and for the normal metal 3 we simply write a transmitted 
electron and hole: 



(31) 



The matching conditions at x = e < give the equa- 
tions: 



Vi(£-)=Vi(e + ), 



(32) 
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We now proceed with the matching condition at x = 
using the theory in Ref£. The single-valuedness of the 
wave function at x = implies that: 



du 



dv 



a + (3 = cu + 
a + (3 = t, 

7 + 6 = cv + 

7 + 6 = t a , 

and the (probability) current conservation implies that 
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du_ + —t , (42) 



iq 



dv_ - -2-t , (43) 
m hf m 



where m denotes the electron's effective mass in metal 3. 
Equations (|34 p - (|3~T|) allow the elimination of the ampli- 
tudes a, (3, 7, S. 

These equations can be applied again to the case of 
normal incidence at E = 0. Defining the Fermi momen- 
tum and velocity in metal 3 as hqp and w 3 = hqp/m 3 , re- 
spectively, we obtain modified results for the amplitudes 
a and b: 
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In the limiting case of a clean junction with the nor- 
mal incidence close to the nodal direction, Z, k — > 0, the 
above expressions simplify to 
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FIG. 1: Scattering amplitudes (|46[) and conductance G = 
1 + \a\ 2 — \b\ 2 for Z, k — as functions of rf for: rj = rf (left 
panel); r\ = 0.575 (right panel), which is the value inferred 
from the effective barrier parameter estimated in Ref.— . 



where r\ = v F /v n , as in equation ([27]) . and 77' = vz/v n . 



tf, 

In order to obtain the perfect conductance for this case, 

we assume no mismatch of light Fermi velocities (77 = 

rf = 1) so that |6| 2 = 0.04 and \a\ 2 = 0.16. The perfect 

conductance in this case is then 



G 



perfect 



1.12. 



(47) 



which should be compared to the well known result 
1 + |a| 2 - \b\ 2 = 2 in BTK theory. The effect of metal 
3 is to "short-circuit" the system by providing a new 
transmission channel for the incoming electron. In this 
case the electron does not need to pick up a second elec- 
tron and penetrate the superconductor as a Cooper pair, 
leaving the Andreev hole behind. Hence the suppres- 
sion of the Andreev reflection process from \a\ 2 = 1 to 
I a. 1 2 = 4/25. The effect of metal 3 is reduced if the ratio 
i>3 jv n decreases and the theory of section [TT| is recovered 
in the limit v^/v n — > 0. In the oposite limit, v^/v n — * 00 
implies a — > and b — > — 1 . Figure Q] shows a plot of 
equations (|4"6"| and the conductance for some choices of 
velocity ratios. 

We note that the result (T4"7| is remarkably close to 
the enhanced subgap conductance observed in CeCoIns, 
where the zero-bias differential conductance^ lies in the 
range 1-1.13. Figure Q] shows that some choices of veloc- 
ity ratios yield conductances in this range. A fit of our 
model to the extensive experimental data on this com- 
pound might settle the model parameters values. 

We also note that the form of equations (|4"4")) - (|4"5"|) shows 
that the conductance cannot be written a sum of two in- 
dependent parallel conductances, one from the heavy liq- 
uid, another from the light one. The interface impedance 
here is a quantum-mechanical effect due to the boundary 
conditions for the wave function at an intersection, as in 
quantum waveguide theory. 
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IV. CONCLUSIONS 

The success of the two-fluid model used in Rcf. 5 
to explain the conductance spectra of the junction 
Au/CeCoIn 5 relies on the assumptions that: (i) only the 
light velocities are important; (it) the density of states 
is decreasing as one approaches the Fermi level; (in) 
the conductance may be obtained as a weighted aver- 
age. Here we have explicitly confirmed the validity of (i). 
Item (ii) is not consistent with the usual approach where 
one considers the lowest band partially filled, but is con- 
sistent if we consider that the lowest band is full and the 
Fermi level is located at the bottom of the higher band. 
In this regime the effective mass is also high and the den- 
sity of states is decreasing with energy. However, only the 
light velocities affect the Andreev reflection. Regarding 
(Hi), we have shown that in general the conductance may 
not be written as a sum since the interface impedance is 
a quantum-mechanical effect which leads to interference 
terms. 
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FIG. 2: Left: the lower partially occupied band E-(\t) of 
a normal HF system resulting from the hybridization be- 
tween the / and c subsystems with energies e = e/ — fi 
and £ = £k — Mi respectively. The c band energy evalu- 
ated at the Fermi momentum, hp, is £f and has a slope 
= (d£/dk)F which is larger than the actual Fermi veloc- 
ity Ti~ 1 (dE-/dk) F = v^e 2 /V 2 . Right: the upper partially 
occupied band E + (k). 
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APPENDIX A: ANALYSIS OF EQ. © 

If A = and the system is less than half-full, all elec- 
trons are in the lower band given by 



£L(k) 



AV 2 



The Fermi level is given by E- = 0, which means that 
0=0 in equations © and (J7J), or: 



f - V2 

?F - 

e 

The Fermi velocity is then given by 



dE_ 
~d~k 



UE_ 



2 hv F 



(Al) 



(A2) 



where Vp = hT 1 (d^/dk) F denotes the c-band velocity 
evaluated at the Fermi momentum Tikp, and is compara- 
ble to that of the light metal. In order for the density of 
states (or mass enhancement) to be high, we must have 
V 2 ^> e 2 . Figure [2] shows the bands in the normal state 
of the HF and the relevant energy scales. 

For finite A, the lowest band of excitation energies 
is obtained by choosing the minus sign in equation J5J, 
which may be written as 



E 2 




(A3) 



where 70 and (3q are the functions 7 and (3 given in ([6])- 
([7]) when A = 0. By expanding in A (assuming V 2 ^> 
e 2 3> A 2 which is valid for instance for CeCoIn 5 where 
the coherence temperature is of the order of 45 K and 
the superconducting critical temperature is of the order 
of 2.3 K) we obtain the gap for excitations at k F '- 



E(k F ) = A(k F ) 



A(k F ) (A4) 



and the excitation spectrum in the vicinity of the Fermi 
level is: 



E(k) 



\ 



hvf (|k| - k F ) 
1 + m 2 



A 2 (k) 



(A5) 



which has the usual form and involves the slow velocity 
of the heavy electrons. 

We may as well consider a situation where the local 
Coulomb repulsion between the electrons is not too large, 
which enables a /-site occupancy nt > 1. This case is 
also represented in Fig. [2l where we consider that the 
Fermi level is located above e/. In this case £f,£ < 0. 
The position of the Fermi level is again obtained from 



(3q = and leads in this case to £f 



-V 



It is 



easy to see that the Fermi velocity is once again of the 
form (|A2[) . as expected of a heavy band. Note however, 
that while for the case when rif < 1 the density of states 
is increasing as we approach the Fermi level, in the case 
when rif > 1 the density of states is a decreasing function 
of the energy in the vicinity of the Fermi level. 
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